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Abstract 

A nonperturbative approach to QCD describing confinement and 
chiral symmetry breaking is discussed. It is based on the path in- 
tegral representation of Green's function of quarks and leads to the 
QCD string theory. The effective actions for mesons and baryons in 
the external uniform static electromagnetic fields are obtained. The 
area law of the Wilson loop integral, the approximation of the Nambu- 
Goto straight-line string, and the asymmetric quark-diquark structure 
of nucleons are used to simplify the problem. The spin-orbit and 
spin-spin interactions of quarks are treated as a perturbation. Using 
the virial theorem we estimate the mean radii of hadrons in terms 
of the string tension and the Airy function zeros. On the basis of 
the perturbation theory in small external electromagnetic fields we 
derive the electromagnetic polarizabilities of nucleons. The electric 
and diamagnetic polarizabilities of a proton are = 10 x 10^^ frrfi, 
pdia _ _g X 10~^ fm^ and for a neutron we find q;„ = 4.2 x 10~^ fm^, 
Pn°' — ~5.4 X 10~^ fm^. Using the A contribution to the paramag- 
netic polarizability of the nucleons, reasonable values of the magnetic 
polarizabilities = (5 ± 3) x 10"'' fm^, = (7.6 ± 3) x IQ-^ fm? 
are estimated. 



1 Introduction 

It is now understood that Quantum Chromodynamics (QCD) (see [1]) re- 
ally describes the strong interactions of hadrons. However, QCD deals with 
nonobserved objects: quarks and gluons. The standard perturbation theory 
in the small QCD coupling constant (a^) uses Feynmann diagrams and is 
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applicable in the ultraviolet region. In the infrared region, which is impor- 
tant for hadron physics, the coupling constant is large, and it is impossible 
to apply perturbation technique. In this case the nonperturbative meth- 
ods should be used. One of the powerful methods is lattice simulations of 
QCD giving data about hadron spectrum, matrix elements, and scattering. 
The second way is to develop analytical nonperturbative methods which al- 
low us to understand deeper nonperturbative phenomena. There are some 
phenomenological approaches for describing hadron physics such as poten- 
tial models [2], method of the effective Lagrangians [3], the QCD sum rule 
approach [4] and others. However, these schemes do not deal with the fun- 
damental QCD Lagrangian. 

One of the important problems of particle physics in the infrared region is 
the confinement of quarks. According to this phenomenon only colourless ob- 
jects (hadrons) are observable. This means that quarks and gluons which are 
present in the QCD Lagrangian are absent in the physical spectrum. Forces 
between quarks and gluons combine them into hadronic colourless states. The 
appearance of such forces with a linear potential is explained by the string 
configurations of gluonic fields named the QCD string. Mathematically this 
phenomenon is expressed by the Wilson area law. Confinement of quarks 
and gluons can be explained by the dual Meissner effect [5] in the frame- 
work of the Abehan- projected method [6]. In accordance with this scenario 
Abelian magnetic monopoles appear after Abelian projecting Yang-Mills glu- 
onic fields. These magnetic monopoles condense to form vacuum with fiux 
tubes (strings) between quarks. The same picture occurs in the Ginzburg- 
Landau theory of superconductivity and its relativistic version - Abelian 
Higgs Model. In these theories the magnetic Abrikosov-Nilsen-Olesen strings 
[7] appear due to the condensation of Cooper pairs (electrons). Thus, chro- 
moelectric strings which ensure confinement in the framework of QCD are 
dual to the Abrikosov-Nilsen-Olesen strings. 

It should be noted that QCD strings are treated as effective and are real 
gluonic flux tubes. This is a distinction of QCD strings in AD which can 
be described by the Nambu-Goto action [8] from the fundamental strings 
[9]. Therefore there is no problem to quantize QCD strings and we avoid 
the problem of conformal anomaly in D 26. Superstring theories were 
introduced to unify strong, electroweak, and gravitational interactions. 

The goal of this review is to outhne the effective string theory which 
appears from QCD in the infrared region, and allow us to make some esti- 
mations of electromagnetic characteristics of hadrons. 
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There is another nonperturbative phenomenon of QCD - chiral symmetry 
breaking (CSB) the QCD Lagrangian with massless quarks. For massless 
quarks, the QCD Lagrangian possesses the chiral U{Nf) x U{Nf) symmetry 
which is spontaneously broken. As a result there is no parity degeneracy 
of hadrons. The real difference of hadron masses with opposite parity is 
few hundred MeV. The consequence of chiral symmetry breaking QCD is 
the presence of the Goldstone bosons in the physical spectrum, which are 
associated with vr-mesons. 

CSB leads to a nonvanishing quark condensate {{qq)) [10]. As a result the 
light quarks {u, d quarks) with current masses — ma — 7 MeV acquire 
the dynamical masses ~ //.^ ~ 320 MeV. This phenomenon is important 
for light pseudoscalar mesons because they possess the Nambu-Goldstone 
nature. Nonzero gluon condensate is responsible for confinement [4]. To 
obtain low masses of pseudoscalar mesons one needs to take into account the 
spin interactions of quarks. 

Some models were proposed for the explanation of CSB. In instanton 
vacuum theory (IVT) [11], CSB was explained by the nonvanishing density 
of quark zero modes. This approach is based on the representation of vacuum 
as an ensemble of classical instanton- antiinstanton configurations. To have 
zero topological charge of the QCD vacuum the number of instantons should 
be equal to the number of antiinstantons. However, a superposition of fields 
with opposite topological charges is instable as instantons and antiinstantons 
can annihilate, and in addition such a superposition of fields is not a solution 
of classical equations of motion. Classical configurations should have high 
density to reproduce the phenomenological gluon condensate [4], and can 
distort the original solutions. Besides, IVT can not explain confinement of 
quarks. 

In [12] CSB was confirmed on the basis of gluon propagator which is the 
solution of Dyson-Schwinger equations. 

CSB was also studied in detail in the framework of Nambu-Jona-Lasinio 
models [13-15] based on four-fermion interactions. 

Recently a new analytical method for studying nonperturbative effects 
of QCD appeared [16-20] called Method of Field Correlators (MFC) which 
is the generalization of the QCD sum rule approach [4]. MFC is based on 
QCD and takes into account both gluonic strings and CSB. The fundamen- 
tal variables of MFC are the gluonic condensate and the correlation length 
of gluonic fields. This method is based on first principles of QCD and is 
supported by experiment and lattice data. It can solve problems which can 
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not be solved in the framework of relativistic quark models. Such problems 
inchidc: (i) the Reggc slopes of meson and baryon trajectories and their in- 
tercepts, (ii) constituent masses of quarks, (iii) radiuses and electromagnetic 
polarizabilities of hadrons, (iiii) high energy scattering, and others. 

In the framework of MFC both non-perturbative effects of strong inter- 
actions (confinement and CSB) can be explained by introducing stochastic 
gluon vacuum fields with the definite fundamental correlators [16,17]. Then 
the linear potential between quarks appears and it provides the confinement 
of quarks. The confinement of quarks does not allow them to be observed, 
i.e. quarks can not move outside of hadrons at large distances relative each 
other. This was confirmed by Monte-Carlo simulations and experiments. 

The Regge trajectories are asymptotically linear in this approach with 
a universal slope [17]. There are different stochastic vacuum configurations 
which are responsible for CSB: instantons (antiinstantons), pieces of (anti- 
) self-dual fields (for example torons or randomly distributed lumps of field) 
and others. The necessary requirement is to have zero fermion modes. The 
condensation of zero modes leads to CSB. 

So the MFC in nonperturbative QCD and the dynamics of zero modes 
give the explication of the double nature of light pseudoscalar particles (pions, 
kaons and others) as Nambu-Goldstone particles and as the quark- antiquark 
system with a confining linear potential. It should be noticed that confine- 
ment prevents the delocalization of zero modes over the whole volume [18] i.e. 
stabilizes the phenomena of CSB. The familiar PCAC (partial conservation 
of axial vector current) theorems and the soft pion technique are reproduced 
in this approach [19]. 

There are some difficulties in evaluating meson characteristics in the gen- 
eral case of the complicated string configuration. Naturally, as a first step, 
we make some approximations and model assumptions to simplify the calcu- 
lations. So here we consider the straight-line string as a simple configuration 
and quarks attached to the ends of the string. Such configurations were also 
studied in [21]. The quark-diquark structure of baryons is employed here. 

Spin degrees of freedom are treated as a perturbation and therefore it 
is questionable to apply this scheme to pions and kaons. For example p 
and K* mesons can be considered here because the energy shift for them 
due to the hyperfine spin interaction is below 100 MeV [2]. Here short- 
range spin-orbital L • S and spin-spin Si • S2 interactions are not taken into 
account. The Coulomb like short-range contribution due to the asymptotic 
freedom of QCD can be easily included. We imply also that in the presence 
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of light quarks the structure of vacuum yields an area law of the Wilson loop 
integral. In the present review we investigate hadrons in external, constant, 
and uniform electromagnetic fields and use the path integral approach. 

It is important to calculate different intrinsic characteristics of hadrons 
on the basis of QCD string theory and to compare them with the experi- 
mental values. It will be the test of this scheme. The mean radiuses (and 
electromagnetic form-factors) and electromagnetic polarizabilities of hadrons 
are fundamental constants which characterize the complex structure of par- 
ticles. These values for some hadrons are known from the experimental data 
and therefore the estimation of them is reasonable. 

The electromagnetic polarizabilities of hadrons a, /3 enter the induced 
electric D = «E and magnetic M = /9H dipole moments, where E, H are 
the strengths of electromagnetic fields. As a result there is a contribution to 
the polarization potential [22] as follows U{a,(3) = -{l/2)aE'^ - {l/2)/3H'^. 
Electromagnetic polarizabilities are fundamental low-energy characteristics 
of strong hadron interactions and therefore they can be calculated in the 
framework of nonperturbative quantum chromodynamics - QCD string the- 
ory 

The review is organized as follows. In Section 2 after describing the 
general background we derive the Green function of quark- antiquark sys- 
tem (meson). The effective action and Hamiltonian for mesons in external 
electromagnetic fields is found in Section 3. In Section 4 we estimate the 
mean size of mesons on the basis of the virial theorem. Section 5 contains 
the derivation of the Green Function of baryons (three-quark system) . The 
effective action for baryons in external electromagnetic fields based on the 
proper time method and Feynman path-integrals is derived in Section 6. 
Section 7 contains the calculation of average distances between quarks in 
nucleons. The electric polarizabilities of nucleons are evaluated in Section 
8. We derive diamagnetic polarizabilities of protons and neutrons using the 
perturbative expansion in the small magnetic fields in Section 9. In section 
10 the method of field correlators and cluster expansion are briefiy described. 
The introduction to the hadron electromagnetic polarizabilities is given by 
Appendix. In the conclusion we made the comparison of our results with 
other approaches. 

Units are chosen such that h — c — l. 
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2 Green's Function of Quark- Antiquark Sys- 
tem 



Here we derive the Green function of quark-antiquark system using the 
Schwinger proper time method and the Feynman path-integrals. Our goal is 
to calculate some electromagnetic characteristics of hadrons. For this purpose 
we need the effective actions for mesons and baryons in external electromag- 
netic fields. The method of the Green functions will be adopted. 

Let us consider the Lorentz and gauge invariant combination of two quark 
colorless system (meson) in Minkowski space [17] 

Xm{x, X, C) = qa{x)rA^ab{x, x)qb{x), (1) 

where qb{x) and qa{x) are quark bispinors, qa{x) = qa{x)'ji; qa{x) is the 
Hermitian-conjugate quark field; 1, 75, 7^, ^Ibln, = -{'i/4:)bn,lu] 

(F^ = Fa, i.e. F^ are Hermitian matrices); 7^ are the Dirac matrices; a, b 
are colour indexes and as usual there is a summation on repeating indexes. 
We imply that quark fields q{x), q{x) possess the definite flavours which form 
mesons. The gauge invariance is guaranteed here by introducing the parallel 
transporter [17]: 

^x,x) = Pexp l^ig J_ A^dz^^ , (2) 

where P is the ordering operator along the contour C of integration, g is the 
coupling constant, = v4^A"; are the gluonic fields; A" are the Gell- 
Mann matrices. Although we have an attractive feature - gauge invariance, 
the function (1) depends on the form of the contour C. As Xm{x, x, C) is 
a gauge invariant object, it obeys the Gauss law on the spacelike surface E. 
The path of integration in Eq. (2) is arbitrary. 

Let the four-points x, x, and |/, be the initial and final positions of 
quark-antiquark, respectively. Two particle quantum Green function is de- 
fined as [17]: 

G{x,x;y,y) = {XM{x,x,C)X^{y,y,C')), (3) 
were Xf^(y, y, C) corresponds to the final state of a meson: 

xtiiy. y, C) = qa{y)^ab{y, y)TAqb{y), (4) 
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where Fa — •JaTaI'I- The brackets (...) mean the path-integrating over glu- 
onic and quark fields: 

{XmX+) = J DqDqDA^ exp {zSqcd} XmX+, (5) 

with the QCD action Sqcd- We imply that the measure DAfj, in the path 
integral (5) includes the well known gauge-fixing and Faddeev-Popov terms 
for gluonic fields [23] . The Minkowski space is used here but it is not difficult 
to go into Euclidean space to have well defined path-integrals. 

Let us introduce the generating functional for Green's function to calcu- 
late the path-integral (5) with respect to quark fields: 

Zl,M-lD,D,e.,{,S,an + .lM,M)+UM}, (6) 

where we introduce the external colour anticommutative sources (Schwinger 
sources) rja, fja- Then the Green function (3) can be written as 



G{x,x;y,y) = j DA^^ 



S 5 



s s 



(7) 



ri=ri=0 

Now it is possible to integrate the path-integral in Eq.(7) over quark fields 
q, q as expression (6) is a Gaussian integral. We may represent the QCD 
action in the form of 

Sqcd = S{A) - J dxq{x) {j^D^ + m) q{x), (8) 

where S{A) is an action for gluonic fields with the included ghost fields, 
Dfx — dfj_ — iQA^^ — igA^] A^^ is the vector potential of an electromag- 
netic field, Q is the charge matrix of quarks, Q — diag{ei,e2, ■■■eNf), e-i 
are charges of quarks, Nf is the number of fiavours; m is the quark mass 
matrix, m = diag{mi,m2, ...mN^) and we imply the summation on colour 
and flavour indexes. Thus we have just introduced the interaction of quarks 
with an electromagnetic field. Inserting Eq. (8) into Eq. (6) and integrating 
with respect to quark fields, we arrive at the expression 

Z[f],r]] = det(-7^D^ - m) exp ^iS{A) + i J dxdyf){x)S{x,y)r){y)^ , (9) 
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where the classical quark Green function S{x, y) is the solution of the equa- 
tion 

{^^D^ + m)S{x,y) = 6(x-y). (10) 

Using Eq. (9) and calculating the variation derivatives in Eq. (7) we find 
the quantum Green function of a meson: 

G{x,x;y,y)^ j det(-7/,D^ - m) exp {^^(A)} 



X tr 



^ASna{y, x)TA^ab{x, x)Sb„^{x, y)^mn{y, 2/) 



-tr 



^mn{y,y)^ASnm{y,y) 



tr 



^A^ab{x,x)Sbaix,x) 



(11) 



The second term in (11) corresponds to the annihilation graphs and it does 
not contribute to the nonsinglet channel. The functional determinant in 
Eq. (9) describes the contribution from the vacuum polarization and gives 
the additional quark loops. In the quenched approximation the fermion de- 
terminant is taken to be unity. Neglecting quark- ant iquark vacuum loops 
(quenched approximation) and omitting the annihilation graph, the Green 
function of a meson (the quark- antiquark system) takes the form (see also 
[17]) 

G{x, x; y, y) = (tr[f a5'„„(7/, x)TA^ac{x, x)Scm{x, y)^mn{y: y)])A: (12) 

were the brackets {■■■)a mean the averaging over the external vacuum gluonic 
fields with the standard measure exp [iS (A)]. 

Let us derive the classical one-quark Green function which is the solution 
of Eq. (10) using the Fock-Schwinger method. Starting with the approach 
[24], the solution to Eq. (10) for the Green function of a quark (in Minkowski 
space) is given by 



S{x,y) = {x I (D + mi)-^ \ y) = {x \ (mi - D){ml - D^)'^ \ y) 



= (?(^)?(y)), (13) 

where Ci and nii arc the charge and current mass of the quark, D ~ 7//-D^, 
-D^ = — ieiA^^ — igAf^, A^ = A^A"; 7^, A" are the Dirac and Gell-Mann 
matrices, respectively; A'jj, A"^ are the electromagnetic and gluonic vector 
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potentials. The inverse operator (ml — D'^^ can be represented in the 
proper time s [24] : 



(ml-D^) ^ = i ds exp {-is [ml - } ■ 



(14) 



Using the properties of the Dirac matrices {7^, ji,} = 2S^i, we find the squared 
operator 

D'^Dl + E^. (eiF^l + gF,J) , (15) 



where 



pei _ o /ie( _ o Ael 



el 



\el 



(16) 



Tj^iy are the spin matrices, and F^l and are the strength of electromag- 
netic and gluonic fields, respectively. Inserting relationship (14) into Eq. (13) 
with the help of Eq. (15) we get 



roo 

S(x, y) — i dslx I (mi — D) 
Jo 



X cxp{-^sK -Dl- ^.uieiF^l + gF,,)]} \ y). (17) 

The exponent in Eq. (17) plays the role of the evolution operator which 
defines the dynamics with the "Hamiltonian" m\ — 

with initial | y) and final {x \ states where s means the proper time. Therefore 
it is convenient to represent the matrix element in Eq. (17) as a path integral 
[25]: 



Six, y)^i dsN / DpDzP (mi - d) exp i , 

Jo Jz{Q)=y ^ ^ I Jo 

- (p^ - eiA'l - gA^ + E^, (e^F^ + gF^^ 



mi 



(18) 



where {t) is the path of a quark with the boundary conditions z{0) — y, 
z{s) — X, D — — eiA^j^ — gAjj^ and P means ordering; A^" is a constant 

which is connected to the measure definition and it will be chosen later. The 
path integration over the momenta can be rewritten in the form (see [17,26] 

N J Dp (mi - exp |i ^ dt p^^z^, - {p^ - CiA''^ - gA^f | 
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N 



I Dp 



exp 



X exp |-i ^ dt (p^ - eiA^j - 
J Dp (mi- ^7m^m) exp |i ^ o?t p^i^ - (p^ - CiA^^ - gA^Y | 
J Dp (mi- ^7^i^) exp |z dt -pi + ^ij + {ciA^^ + ^yl^) i 



(19) 

In Eq. (19) we used the integration by parts (see [26]) and made a continuity 
of shifts P;ti — > + eiA^J; + gA^ and then p^ ^ Pf^ + The constant N 

in Eq. (19) is defined by the relation 



TV J Dpexpl^^-zJ^dt{pl)^ 



(20) 



Taking into account Eqs. (19), (20) we find from Eq. (18) the Green function 
of the quark: 

roo /■z{s)=x ( % \ { 

S{x,y) = i ds Dz [mi - --fi^z^{t)) Pexp \i dt 

Jo Jz(0]=y \ 2 J y Jo 



z{0)=y \ 2 

+ (eiA^l + gA^ z, {t) + S^, (eiF^^ + gF, 



(21) 



3 Effective Action for Mesons 

To get the effective action for mesons in the external electromagnetic fields 
we use the Green function (21) of the quark possessing spin in the Minkowski 
space which can be represented as 

roo rz(s)=x / n \ r /■« r 1 

S{x,y)^i ds Dz (mi - ^ ^^^^ ^ — ' -^..n ^2 

^0 Jz(0)=y \ 



^^t^^t^ (t) jPexpii -z (t) - mi 
2{o)=y V 2 / I io L4 



+ eiZ,{t)A^;iz) + E^. (eiF^l + gF,,) 



dt\^{x,y), 



(22) 



where $(a;, y) is the path ordered product (2). Inserting Eq. (22) for a quark 
{S{x,y)) and antiquark {S{y,x)) Green functions into Eq. (12) we find the 
expression for the Green function of a meson: 



rco roo ^ rz(s)=y _ / ^ 

G{x, x; y, y) = tr ds J - ^^^^ ~ 2^1^^/^ (^) 
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x(Ps-PAexp<^ic/E 



F^u [z] dt 



Va exp \ -igT,^^ / F^^ {z) dt 



]W{C))^ 



X 



z{s)='X 

m=y 



Dzim-z- -^t,z^{t) 



X expji -zlit) -ml + eiz^{t)A''J;{z) + ciE^^F^l {z) 



dt 



+ i 



-zlit) -ml + e,z,{t}A^^{z) + e.E^^F^l (z) 



dt 



(23) 



where Ps and Pa are the ordering operators of the spin matrices E^,^ and 
gluonic fields, respectively; ei, 62 are charges and mi, m2 are current masses 
of the quark and antiquark; Zi^{t), z^if) are the paths of the quark and 
antiquark with the boundary conditions -2^(0) = x^, — Un-i ^m(O) — Vn^ 

and z^(t) = dz^/dt. As compared with [17, 27] we added the 

The 



interaction of charged quarks with the external electromagnetic fields, 
gauge - and Lorenz - invariant Wilson loop operator is given by 

1 



W{C) 



N< 



-trP exp <^ ig 



c 



c 



(24) 



where Nc is the color number, and C is the closed contour of lines xx and 
yy connected by paths z{t), z{t) of the quark and antiquark. The Wilson op- 
erator (24) contains both the perturbative and nonperturbative interactions 
between quarks via gluonic fields with the QCD coupling constant g, and 
it is the amplitude of the process of creation and annihilation of quarks and 
antiquarks. 

Using non-Abelian Stokes theorem 

as 



1 



W{C) = ^irPexp 



ig 



, 17], the Wilson loop is rewritten 

(25) 



F^^,{z,zo)da,nAz 



where E is an arbitrary surface bounded by the contour C, and Fi^,^{z, Zq) 
given by 



IS 



F^u{z, zo) = ^zo, z)F^^{z)^z, zo). 



(26) 



The coordinate zq and path of the transporter ^{z, zq) belongs to the surface 
E, and they are arbitrary [17]. The contour C can be parametrized by the 
four vector z^{C,), where two dimensional coordinate ^ = (^1,^2) induces the 
metric tensor of the surface E: g°^l{C} = [daZ^{C}\[df,Zi,{!^)] {da = d/d^a)- The 



infinitesimal element of the surface E is given by 



d(Tf,u{z) 



ah ab 



(27) 
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where e"^ = —e^"',e^^ — 1. In accordance with the approach in [17], spin 

interactions can be treated as perturbations, e.g. for p and K* mesons but 
not for Nambu-Goldstone it- and iT-mesons. To construct the expressions in 
spin interactions we write the relationship [17] 



(exp \ igT, 



f F^, {z) dA\w{C))A 



exp < E 



Jo 



dt- 



where Sa/^i, (t) is the surface around the point (t). The zeroth order in spin- 
orbit and spin-spin interactions corresponds to neglecting the terms ejE^,^F^^, 
in (23). 

To receive the effective action for mesons from Eq. (23) one needs to esti- 
mate the Wilson loop integral (24) which carries very important information 
about an interaction of quarks. 

Let us consider the case when the distance between quarks is greater than 
the time fluctuations of the gluonic fields (r ^ Tg). The Monte-Carlo calcu- 
lations [29] gave Tg ~ 0.2 -^ O.Sfm. We imply that the characteristic quark 
relative distance is r ~ 1 fm. This assumption is confirmed by the lattice 
data [30] and by the calculation of the quark-antiquark relative coordinate 
[31]. 

The average Wilson integral (24) at large distances in accordance with 
the area law can be represented in the Minkowski space as 

{W{C))A = exp{iaSmin): (29) 

where a is the string tension and Smin is the area of the minimal surface 
inside of the contour C. This behavior of the Wilson loop integral is the 
consequence of the nonperturbative interactions of a quark and antiquark, 
and corresponds to the hnear potential. The area law of the Wilson loop 
is the criterion of confinement [32]. In this case the gluonic field between a 
quark and antiquark forms a string, and the string tension a is the energy 
of a string per unit length. In accordance with the lattice data (see e.g. 
[17]) a ~ 0.2 GeV^. It is impossible to calculate a dimensional parameter a 
using the perturbative expansion in the small coupling constant g. There is a 
dimensional parameter of QCD in the high-energy limit - the scale parameter 
Aqcd ^ 200 MeV: 
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where // is a normalization mass. According to the phenomenon of dimen- 
sional transmutation, dimensional parameters are proportional to some power 
of Aqcd- It is easy to see from Eq. (30) that all coefficients of the expansion 
of Aqcd in powers of are zero. Therefore, a is proportional to Aqqj^ and 
the QCD strings are described by the nonperturbative approach. 

In the case of perturbative interactions at small distances (high energy), 
the Coulomb potential dominates, which docs not confine quarks. QCD is 
an asymptotically free theory and at high energy the coupling constant is a 
small parameter. In this regime the gauge field is not concentrated in the 
tube (string) but fills the whole space. Then the asymptotic form of the 
Wilson loop (for a smooth contour C) for this case is given by the perimeter 
law: 

{W{C))a = exp{t^L{C)}{WiC)ren)A, (31) 

where the appearance of the finite vahic {W{C)ren) is due to the renormal- 
ization; 7 is a dimensional constant, and L{C) is the length of the contour 
(L(C) = Jq ds^z'^{s), Zn{0) = -2/i(l)). When the distance between quarks 
r < 0.25 /m, the system is in the Coulomb phase and we have the perime- 
ter law behavior of the Wilson averaged integral (31). The confining phase 
(nonperturbative) occurs at large distances and gives the area law (29). For 
wide class of contours without self-intersections, the Wilson loop integral is 
multiplicatively renormalizable [33] and is given by 

{W{C))a = {W{C))pert{W{C))nonpert, (32) 

where the first factor in right-hand side of Eq. (32) describes the pertur- 
bative contribution in the Wilson integral (31), and the second describes 
confinement. 

The surface Smin in Eq- (29) can be parametrized by the Nambu-Goto 
form [8] 

Srrnn = dv dP^{w^w'^Y - W^Ul'l, (33) 

•J 

where = dw^j^/dr, w'^ = dw^^/dp. 

Using the approximation [17] that the coordinates of the string world 
surface Wh{t,(3) can be taken by straight lines for the minimal surface we 
write 

w,{t, P) = z,{t)P + z,{t){1 - P), (34) 

where r is implied to be the proper time parameter for both trajectories 
T = {tT)/s — {tT)/s. Thus we ignore vibrations of the string. 
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For uniform static external electromagnetic fields we have the represen- 
tation of the vector potential through the strength tensor F^l: 



K'i^) = If^^,, = IK'',- (35) 

The paths z^, are expressed via the center of mass coordinate Rfj, and the 
relative coordinate r^u [17], 

s + s s + s 

with the boundary conditions for i?^(T), r^{T): 

-K^iiUj = ■ , Ki^[I ) = ■ , 

r^{0) =yt^- Vi,, r^{T) =x^- x^. (37) 

The integration with respect to z^, z^ in (23) is replaced by the integration 
over new variables R^, r^. As r is a common time for the quark and antiquark 
(the time of the meson) the parametrization = (z,ir), z^ = (z,ir) is 
possible [17]. This leads to the constraints: Ro{t) — r, ro(r) = 0. In 
accordance with the approach [17] we introduce the dynamical masses /ii, 112 
by relationships 

/^i = ^, = |. (38) 

Replacing the integration with respect to s, s in (23) by the integration over 
d^i and dfi2 with the help of Eqs. (29), (33) - (36) we find [31] the two-point 
function in the zeroth order in spin interactions 

G{x, X- y) = £ DRDr exp {tS^ff} (39) 

with the eflFective action 

'^^^^i '^n-2^-2;^+2^^^ + ^^)^^/^^ 



J 



(40) 
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where u^^ = — /ii/(//i + /i2)]rju, /i = A*iA*2/(a*i+A''2) is the reduced mass 

of the quark- ant iquark system, e = ei + 62, g = Ci — 62. As a first step we are 
interested here in the spinless part and therefore the preexponential terms 
(nil — ^Ifi^fi (j^2 — ^Ifj. {t}^ and the constant matrices Tj^^F^jI were 
omitted in Eq. (39). The expression (40) defines the effective Lagrangian for 
mesons in external uniform static electromagnetic fields in accordance with 
the formula Seff = Jq Cej/dr. The expression (40) looks like nonrelativistic 
one at F^^, = 0, but it is not. The author of [17] showed that the relativism is 
contained here due to the fl dependence and the spectrum is similar to that 
of the relativistic quark model. 

The mass of the lowest states can be found on the basis of the relationship 

[34] 

J DRDre^p{iSeff} = 
A*i + A*2 A*i + A*2 

(41) 

where the mass Ai{fii, ^2) is the eigenfunction of the Hamiltonian. After that 
the Green function (39 ) is derived by integrating (40) over the dynamical 
masses /^i, //2- In accordance with [17] we estimate the last integration on 
d/ii, d^2 using the steepest descent method which gives a good accuracy when 
the Minkowski time T ^ 00. To have the correct formulas, it is necessary to 
go into Euclidean space and return into Minkowski space on completing the 
functional integration. We use this procedure. 

The last term in (40) can be represented by the relation 

12 



J 



a j dj3^ 



L 60 + 261 1/?-^^ +62 , (42) 

\ \ III +1^2 J V 1^1+1^2/ 

where 

60 = (r^^i?^) - rlRl, 61 = (r^r^) (r^R^) - rl (f^Ry) , 62 = (r^r^)^ - rjr^. 

In the pure potential regime at low orbital excitations of the string when 
the orbital quantum number / is small, expression (42) is equal to o\fx^ . As 
the equalities Rq{t^ — r, ro(T) = are valid, the dynamical quantities are 
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3-dimensional. Prom Eqs. (40), (42) using the standard procedure we find 
the canonical three-momenta corresponding to the center of mass coordinate 



Rn and the relative coordinate r^: 



—a 



dp 



dC 



{r-i^Ri,) r-k - rlRk + {p - ^) [(r^r^) - r^r^ 



—a 







dp^ 



dr'k 



Ml 



/i4 + ^Klr. + -/iR, 



) 



Tk - rlRk 



+ 



2 r 



bo+2b,{p-^)+h{p-^y 



(43) 

The Hamiltonian 7i = TTkTk + ^kRk — J^eff found from (40) with the help of 
(42), (43) takes the form 



^m^^ ^ /^j^2 + ^.2 _ e^ER) - l(Er) 

2 



2/^1 2/^2 

+ (T / dp 
Jo 



6o + 2h (p - 



+ b2 



IJ.1+IJ.2) 



(44) 



so that the equation for the eigenvalues is given by 

n^^ M{^Xl,^x2)^. 



(45) 



As the canonical momentum 11 corresponding to the center of mass coordi- 
nate is a constant due to the conservation law, i.e. 11 = const, we can put 
Rk — 0, R4 — i. At this choice parameters occurring Eq. (44) become (at 
r4 = 0) bo = r^, bi = 0, 62 = — (r x i")^- Then the last term in Eq. (44) giving 
the contribution to the energy due to the string is given by 



n 



string 



■ C dp 
Jo 



r^-(rxr)^(/3-^y 
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The expression Hstring makes sense here when the value under the root 
squared is positive. This condition is rcahzed at (r x r)^ < 2r'^ (if fii = /J2), 
i.e. for low orbital quantum numbers /. After the integration over the pa- 
rameter P we get 



string 



arcsm 



Hi / (r X r)2 



A*i + V r 



,2 



+ arcsin 



fi2 /(r X i")^ 



(46) 



In the potential regime at low orbital quantum number when (r x r)^ -C 
Eq. (46) is replaced by 

"^string = CT^Vr^ 

and it describes the linear potential which guarantees confinement of quarks. 

The terms contained the strength of the electric field in Eq. (44) describe 
the interaction of the dipole electric moment d with an external electric field. 
Using the definitions we have 

^(ER) + |(Er) = l(eiri + e2r2)E = dE (47) 

and the interaction energy of the electric dipole moment with a uniform static 
electric field is U = — dE. 

In the center of mass system when R = const and at / = 0, E = 0, H = 0, 
taking into account Eq. (43), we find (see [17]) from Eq. (44) the expression 

The Hamiltonian (48) describes two quarks which are connected by the non- 
rotating string. Finding extremum of Hq in mi and ■m2' dHo/d/ii — 0, 
d7io/Oii2 = 0, one arrives at 



Hi = ^p2 + m? + + ml + a\f^. (49) 

Thus we come to the Hamiltonian of the relativistic quark model (RQM) (see 
[35]). In the approach considered /xi and ^2 are dynamical masses and mi, 
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777.2 a-re current masses of quarks. In RQM Eq. (49) is used for dynamical 
masses mi, m2 of the order of 200 MeV and arbitrary orbital number /. So 
the Hamiltonian (49) ignores the rotation of the string and gives the slope 
of the Regge trajectories l/Scx [17]. Besides, in RQM to get the mass of 
p-meson one needs to subtract 700 800 MeV. Such large negative constant 
introduced by hand in RQM has nonperturbative nature. The appearance 
of this constant is due to selfenergies of quarks and can be explained in the 
framework of CSB in QCD string approach [17]. 

At large orbital quantum number I the Hamiltonian becomes [17] 

nj = 27rayjl{l + l), 
so that the Regge slope is equal to 1/2'na. 



4 Mean size of Mesons 

Now we find the wave function of the ground state and estimate the mean size 
of mesons. The equation for the eigenfunction $o of the auxihary "Hamilto- 
nian" 

~ ml ml (/xi + II2) 

rCo^ no - 

2/xi 2/X2 2 

is given by 

(^p' + aoV^^o = e(/x)$o, (50) 

so that Tio^o = e(A*)$0; where e(/i) is the eigenvalue. We can apply equation 
(50) to the leading trajectories with light quarks with masses mi = m2 = m, 
/xi = yU2 = /U (/i = /i/2) for p mesons and when mi 7^ m2, fJ^i ^ fJ^2 for K* 
mesons [36] . In quantum theory instead of the path integration in r we can 
use the replacement pk — > —id/drk- Then Eq. (50) becomes 



1 92 



2/1 drl 



+ (7oV^)$o = e(/^)^o, (51) 



Equation (51) gives the discrete values of the energy e(yu) due to the shape 
of the potential energy. The numerical solution of Eq. (51) was obtained 
in [40]. It is useful to find the solution to equation (51) for the ground 
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state in analytical form. After introducing the variables pk — (2/i(To)^/^rfe, 
e(/i) = {2il)-^l^al'^a{n) [17], Eq. (51) becomes 

^-^ + pVo(p)=a(n)$o(p), (52) 



V 9p\ 

where p = ^p? + pi + Pi- The solution to Eq. (52) may be chosen in the 
form $o(p) = R{p)^im{d 1 'P) 1 where Y/m(6', 0) are spherical functions. After 
setting the variable R{p) = x{p)/P '^^ come to the equation for the radial 
function 

X\p) + (a{n) - p - ^i^) x(p) = 0, (53) 

where x"(p) = (9^x(p)/'9p^, and / is an orbital quantum number. The 
solutions to Eq. (53) for the ground state / = are the Airy functions 
Ai{p — a{n)), Bi{p — a(n)) [41]. The finite solution to Eq. (53) at p — >• oo 
(/ = 0) is 

x{p) = NAi{p-ain)). (54) 
The constant N can be found from the normalization condition 

X\p)dp^l. (55) 



The requirement that this solution satisfies the condition x(0) = NAi{—a{n)) ■ 
gives the Airy function zeroes [41] a(l) = ai = 2.3381, a(2) = 02 = 4.0879 
and so on. The main quantum number n — + Z + 1, where rir is the 
radial quantum number which defines the number of zeroes of the function 
x(p) at p > 0. For the ground state we should take the solution (54) at 
a{n) = ai (here = 0, / = 0, n = 1): 

Xoip)^NoAi{p-ai). (56) 

For the excited states it is necessary to choose the corresponding value of 
n — rir + I + 1- 

Now let us estimate the mean-squared radius for the state which is de- 
scribed by the function $0 (the solution of Eq. (50)). Multiplying Eq. (50) 
by the conjugated function $q and integrating over the volume we find the 
relations 

(T) + {U) = 



/ 

Jo 
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(T) = / ^dl^dV, {U) = ao 1 V^^^dV. (57) 

It is seen from Eqs. (57) that the mean potential energy (U) — ao (^Vr^'j is 

connected with the mean diameter (^Vr^^ (because r is the relative coordi- 
nate and quarks move around their center of mass), which defines the size of 
mesons. In accordance with the virial theorem [42] we have the connection 
of the mean kinetic energy with the mean potential energy: 

2{T)^k{U), (58) 

where k is defined from the equality U{\r) = \^U{r). In our case of the 
hnear potential k—1 and from Eqs. (57), (58) we get 

{U)^\e{il)^\{2il)-y'ay\{n). (59) 

The use of the steepest descent method for the estimation of the integration 
m. jj. (at H = 0) leads to the conditions [17]: 

where the mass of the ground state A4{iii, 112) is given by [see Eqs. (44), 
(45)] 

Here we consider the more general case as compared with [17] when /xi 7^ 
A*2 {mi 7^ 777-2). This case is reahzed for K* mesons. It is assumed that the 
current mass of u,d-quarks (m^ = 5.6 ±1.1 MeV, rud = 9.9 ±1.1 MeV [43]), 
mi is much less than the dynamical mass {^i ~ 330 MeV), i.e. mi <^ /ii 
and the mass of s-quark m2 (m^ = 199 ± 33 MeV [43]) is comparable with 
III but m2 < /^i. Using these assumptions we neglect the term m\/2iii in 
Eq. (61) and from Eqs. (60) have the equations 

(2/i(To)^/^a(n) = 3/x?, 3m^ ± (2/i(To)^/^a(n) = 3/x^. (62) 

Prom Eqs. (62) we arrive to the expression for the dynamical mass 1x2 (for 
s-quark) : 

A*2 = \/lA + ml. (63) 
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To find /Xi the perturbation in the parameter m"^/ fJ'f will be assumed. Using 
the relation fi2 — + "^^2/(2/^1)) which is obtained from Eq. (63) and the 
definition of the reduced mass jl = Aii/^2/(A*i + A*2) from Eqs. (62) we arrive 
at the equation 



3/4 / 2 
' 1712 



a[n) \ I ^ 1712 



/^i^V^(^) (1 + ^ • (64) 



In zeroth order we come to the value /xq = A'-i^'' = \/^{ci{n)/3)^^^ [17]. The 
next order gives the relationship 

In a particular case m2 = we arrive at Hi — ^2 — pi'O — v^(^('^)/^)^^^ 
[17]. The value of the string tension ctq = 0.15 GeV'^ was found from a 
comparison of the experimental slope of the linear Regge trajectories, a' = 
0.85 GeV^^, and the variable a' = 1/8(To [17]. It leads for the lowest state 
rir ^ 0,1 ^ 0,a(l) = 2.3381 to the value //q = 321 MeV [17]. This means 
that for p-mesons when mi = m^, 1712 — rud we have the dynamical masses 
of M, (i-quarks /Xi = /i2 = Ato- Foi' -ft'*-mesons using Eq. (63) and m2 — rUs — 
200 MeV [43] from Eq. (65) we get the reasonable values 

^ll ~ 337 MeV, //2 ^ 392 MeV, fi ~ 181 MeV. (66) 



Inserting the equation {U) = {aovr'^j into the left-hand side of Eq. (59) 
produces the expression 

V^) = ^(2/i<7o)-^/^a(n). (67) 
From Eqs. (63), (65) using the first order in the parameter m\/ ii\ we find 

2A.^„(l + g) f^„ = V^f^My'"|. (68) 

Eq. (67) with the help of Eq. (68) gives the approximate relation for the 
mean relative coordinate 



a[n] \ 



0-0 



1 + 



3m2 / 3 \'/'' 



8crn V airi] 



-1/3 

(69) 
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For p-meson putting m2 = in Eq. (69) we arrive at 

3/4 



rn = ^m) . (70) 



The same expression (70) was found in [31] using another method. With 
the help of the definition of the center of mass coordinate we can write the 
approximate relation for the mean charge radius of p-mesons 0: 



r2) ^ 1 (v^) . (71) 
At (To = 0.15 Gel/2 ^^7^44] ^(^^^ ^ 2.3381 Eqs. (70), (71) give 



rj) ~ 0.42 fm ( (V^) = 0.84 fm) . (72) 



The value (72) characterizes the radius of the sphere where the wave function 
of the p meson is concentrated (remember that r is the distance between 
quarks). We know only the experimental data for vr^ mesons which have the 
same quark structure as mesons [45]: 

For calculating the relative coordinate of quarks for K* mesons we should 
use Eq. (67) or (69) with the conditions pi = pu and p2 = Ps (Eqs. (66)). 
As a result formula (67) gives the value of the mean relative coordinate of 
K* mesons: 

v^) = 0.79 fm. (73) 



Ik 

With the help of this relation we can estimate the mean charge radius of K 
mesons 



'rl^ ~ (v^\ = 0.54 (v^\ = 0.43 fm. (74) 

Pl+ P2 ^ ' ^* ^ ' ^* 

The experimental data of the mean charge radius of mesons having the 
analogous quark structure as K* mesons are 



(r|±) = (0.53 ±0.05) /m [46], 



The relationship ^ (r^) ~ ^ v ) is confirmed by the numerical calculations 
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(r^± ) = (0.34 ±0.05) fm^ [45] , 
and for neutral mesons [46] 

7(4^= (0.28 ±0.09) fm. 

Expression (69) gives a value for the charge radius of ir*-mesons, that is 
close to experimental value for /(''''-mesons; experimental data for the charge 
radius of A'*-mesons is not available. 

The first perturbativc one-gluon exchange contribution to the Hamilto- 
nian determines the spin-spin correction such as the Breit-Fermi hyperfine in- 
teraction [19]. The spin-spin interaction is important to explain the Nambu- 
Goldstone phenomenon which takes place for I — s — channel. To take 
into account spin-spin and spin-orbit interactions one needs to use the gen- 
eral expression (23) and calculate averaged Wilson integral with insertions 
[17]. 

Using the Hamiltonian (44) we estimated in [31, 36] the electromagnetic 
polarizabilitics of mesons on the basis of the expansion in small electromag- 
netic fields E, H. The same procedure will be used bellow for nucleons. 

5 Green's Function of Three-Quark System 

Now we consider baryons as a thrcc-quark system. Let us consider the 
Lorentz-covariant and gauge invariant combination of a three quark, colorless 
system (baryon) [17] 

Xb{x, y, z, Q) = eabc [$(^o, x)q{x)l [$(Zo, y)q{y)], [$(Zo, z)q{z)l , (75) 
where q{x) is a quark bispinor; a, b, c are colour indexes so that 

mZo,x)qix)l^^^AZo,x)qA^) 

and Cabc is the Levi-Civita symbol (ei23 = 1). The parallel transporter 
^{Zo,x) is given by Eq. (2). The points x, y, z are the initial coordinates 

of three quarks with fields q{x), q{y) and q{z), and the point Zq is arbitrary. 
Later, the position of Zq will be defined by requiring to have the minimal 
area for the world surface of three-quark system [17]. 

The two particle quantum Green function of a baryon is given by [17]: 

G{xyz, x'y'z') = {Xb{x, y, z, Q)XBix', y', z', CD), (76) 
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were x', y' are final positions of three quarks; Xb{x' , y', z\ C^) corresponds 
to the final state of a baryon: 

Xb{x', y', z', CI) = emnk Wmx\ Z',)\^ W)^{y\ Z'^]^ 4)]. , 

(77) 

where q = q'^j^. Using the generating functional (6) the Green function (76) 
takes the form 



G(xyz,x'y'z')^ J DA^ 



)^bb'iZo,y)^cc'{Zo,z) 



X^rn'mix', Z^)$n'n(y', Z'o)^k'kiz' , Z'^) 

X S^/ Sf]a' {x)Sf]h' {y)Sf]c' {z)5ri^> {x')6r]n' (|/')%' {z')Z[f], t]] 



T]=T]=0 



(78) 



In Eq. (78) Zq and Z'q are the initial and final positions of the string junction, 
respectively. The total surface which consists of world motions of quarks and 
the path of the string junction must be minimal. This requirement defines 
the path of the string junction [17]. Because the path-integral in Eq. (78) is 
a Gaussian, it is integrated over quark fields q, q . 

Taking into account Eq. (9) and calculating the variation derivatives in 
Eq. (78) we find the quantum Green function of a baryon: 



G{xyz, x'y'z' 



DA^ det{- J f,Df, - m) exp {iS{A)] eabc^mnk 



X Sbm{y-i^)San{^-iy)^ck{Z)Z) S'^mi^ i ^ ) ^bniV ) V ) ^cki^ i ^ ) 

+ St^{z, x')Sl{y, y')S:,{x, z') - St^{z, x')Sl{x, y')SUy, z') (79) 

'^^ami^i ^')Sfni^^ y')Sbkiy^ ^ ) ~ ^brnhji ^ )Sfni^^ y')Saki^^ 

where we introduce the following notation for the covariant Green function 

S^JX,X') = $„„,(Zo,x)5„,„.(^,^')^m'm(^',4)- (80) 

As in the case of mesons the functional determinant in Eq. (79) gives the 
contribution of the additional quark loops. The presence of different terms 
in Eq. (80) is connected with the permutations of quark fields because the 
quantum Green function being considered. As different terms in Eq. (79) 
have the same structure, we consider in detail only one term. Neglecting the 
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functional determinant we find tfie approximate expression for the baryon 
Green function 

Gi{xyz,x'y'z') ^- j DA^^ex.^ {iS{A)} eabc^mnkStm{^^^')Stn{y^y')S%{z, z'). 

(81) 

Expression (81) is the basic formula for deriving effective action for baryons. 



6 Effective Action for Baryons 

We consider baryons in external electromagnetic fields. Inserting Eq. (21) 
into Eq. (80) we find 



oo rz{s)=x / ^ 



X exp<j I I dt 



\iUt)-mi + e,z,{t)A^^{z) 



(82) 



where the contour in Eq. (82) consists of lines between Zq, x and Z'q, x' 
and path Zi^{t). Using the expression (82) for each quark, the baryon Green 
function (81) becomes 

Gi {xyz, x'y'z') = i II ds, jUDz"'' (m, - ^7„i<-'> ((, 



xp.nexp{/; 



X exp{^ E £ dt, {z^\t,)f - m] + e,zf{t,)A^l{z^^) 



(83) 



where j = 1,2,3; ej is the charge of the j-th quark; the boundary condi- 
tions 4i)(0) = a;;, z^^^s,) = x, , z(^\0) = y'^, zf\s,) = y,, z(^\0) = 
zj^\s3) — z^ are used here and the Wilson loop is given by (see [17]) 

{W{C,CyC,))A = eabcemnkihcAZo: ^o)) 

\ / am 
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X {^cAZo,Z',))J^c.{Zo,Z',))Ja. (84) 

We took into account relations such as (28) (see [17]). Contours Cj;, Cy, Cz 
correspond to three quarks which have paths z^^\ z^^ and masses mi, 
7712, fns, respectively. Relationship (83) is the generalization of one [17] for 
the case of quarks placed in external electromagnetic fields which possess 
spins. 

We imply further that the average distance between quarks (r) is greater 
than the time fluctuations (in units c = ^ = 1) of the gluonic fields Tg-. 
(r) > Tg. This condition is valid not only for asymptotic baryon states of 
the Regge trajectories with large angular momenta of the baryon but also 
for lower baryon states. The asymptotic of the average Wilson loop integral 
obeys then the area law and is given by (in the Minkowski space): 

{W{C,CyCz))A = exp {-la {Si + S2 + Ss)} , (85) 

where a is the string tension and Sj (j = 1,2,3) is the minimal surface 
bounded by the trajectories of the quarks qj and string junction Zq. The path 
of the string junction is defined by the requirement that the sum 6*1 + 5*2 + 5*3 
is minimal [17]. Following [17], new variables are introduced: 



hT t2T hT T 



Si S2 S3 — 2s,' 



where r means the proper time for every quark and yUj (j = 1,2,3 ) is the 
dynamical mass of the j— th quark. Using Eqs. (85), (86) from Eq. (83) we 
arrive at 



G,{xyz, x'y'z') --i^ Tli^ fli Dz^^ {rn, - i/.,7.i?(r)) 



I V pel 



-ta{Si + S2 + Ss)y (87) 



The integral in the last exponential factor of Eq. (87) represents the effec- 
tive action for three quark system (baryon) taking into account the spins of 



26 



quarks. So T is the time of the observation, r is the proper time of quarks; rrij 
and the current and dynamical masses of j— th quark. It follows from 

Eq. (87) that there is integration over dynamical masses /ij to get the Green 
function of the three quark system. Pre-exponential factors in Eq. (87) al- 
low us to calculate in principle the spin-spin and spin-orbital contributions 
to the effective action. As a first approximation (see [17]) we neglect the 
short-range spin corrections and consider therefore scalar quarks. The terms 
ejTi^i^F^l which describe the interaction of the magnetic field with the spins 
of quarks will be omitted. With this assumption we arrive at the effective 
action for baryons in external electromagnetic fields 



B 



+ e.-i?'Hr)^?(^(^')) 



a {Si + S2 + S3) . (88) 



The case when electromagnetic fields are absent was considered in [17]. The 
terms which describe the interaction of quarks with electromagnetic fields 
are essential for us because we are going to calculate electromagnetic charac- 
teristics of nucleons. It is convenient to introduce new variables i?^, and 
77^ instead of zj^^ in accordance with relationships [17]: 

/ \ 1/2 / \ 1/2 

where i?^ is the center of mass coordinate of a baryon; and 7]^ are relative 
coordinates of quarks, M — /ii + 112 + IJ'3 is the sum of dynamical masses 
of quarks. The arbitrary mass parameter fi in Eq. (89) defines the scale 
of relative coordinates and rj^. After the substitution (89), the measure 
YljDz^^^ transforms into DRDrjD^ in path- integral (87). 

Let us consider the uniform and constant external electromagnetic fields. 
Then the vector-potential of electromagnetic fields can be represented as 
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Inserting Eqs. (89) into Eq. (88) we find the effective action for the three 
quark system (see [37]) in the form 



B 



= / 

■/o 



dr 



■(7(>5i + >S2 + -53)+AS, (90) 



where we introduce parameters: 



(91) 




A 



M 



A*i + A*2 



(ei + 62) /X3 _^ 63 (/ii + //2) 



A*! + A*2 



(92) 



and e = ei + 62 + 63 is the charge of a baryon. As a particular case, when 
electromagnetic fields are absent {AB = 0) we arrive at the action derived 
in [17]. It follows from Eq. (90) that the center of mass coordinate R^ is 
separated from relative coordinates and rj^ and /j, plays the role of the 
mass of the r/^ excitations. Following [17], the straight line approxima- 
tion for strings and the asymmetric quark-diquark structure of baryons will 
be assumed. The asymmetric configuration (sec also [35,38]) means that two 
quarks q^^^ and g*^^-* are near each other and quark q^^^ is farther from them. 
This case is preferable [35,38] and the slope of linear baryon Regge trajecto- 
ries is the same as for mesons [17]. Then » ^/rf (^^ = ^1+^2 + C3) ^nd 
the coordinate 1]^ can be ignored. We neglect therefore the surfaces Si, S2 
and assume for 5*3 the following expression [17]: 



Ss = b[ dT^\ 

Jo ^ 



(93) 
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where 

Equation (93) takes into account the confinement of quarks and gives the 
hnear potential between quarks. Using the definition B = Jq drC, where C 
is the Lagrangian, from Eq. (90) we arrive at the effective Lagrangian for 
baryons 

^eff = f + f a - E ^ - -b^' + (95) 
Here we neglect the coordinate r]fj_ and introduce the notation: 

eR^R, + Xi^c^u + 7 (^mC. + ii.Ru)] , (96) 



Lagrangian describes the electromagnetic interaction of the string. It 
follows from Eq. (89) that 



As the second term in Eq. (97) is small, the coordinate is proportional to 

the "distance" between quark g*^'^^ and the center of mass of quarks q'^^^ and 
which form a diquark. At the large time T limit ^4 = 0, i?4 = it [17] 
and Lagrangian (96) takes the form 

= e (RE) + 7 (^E) - ]^tmnkHk (eRmK + Kmin + liRrain + LRu)) , 

(98) 

where the electric field — iF^^ and magnetic field — {l/2)ekmnFmn- 
To clarify the physical meaning of the terms in Eq. (98), let us consider the 
dipole moment of quarks. Using the definition of the electric dipole moment 
and Eq. (89) we have 

j 

So, first two terms in Eq. (98) ( neglecting the coordinate rj^ ) describe 
the interaction of the dipole moment of quarks with the electric field in 
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accordance with the expression for the potential energy: U — — (dE). The 
magnetic moment is given by 



-1 ^ 



p . M) f (i) 



2 ^mnk 



+ PVmVn + S [RmVn + VmRn) + 



(100) 



It follows from Eqs. (91), (98) that there is an interaction of the magnetic 
field with the magnetic moment in such a way that the interaction energy is 
U — — (mH). So Lagrangian (98) describes the interaction of electric and 
magnetic moments of baryons with electric and magnetic fields, respectively. 



7 Mean Relative Coordinates of quarks in 
Nucleons 

Now we estimate the mean size of baryons. At the large time T limit ^q= 
(■^4 = "^Co)) -R4 = [17] and therefore only three dimensional quantities Rk 
and are important. From Eq. (98) wc find three momenta corresponding 
to the center of mass coordinate Rk and relative coordinate '■ 

dC ff ■ 1 

Hfe = — ^ = MRk - -eknmHm {eRn + '^lin) , 

oRk 2 

dCeff ■ 1 

TTfe = —4^ = H^k - 7,^knmHm (A^n + '2'yRn) . (101) 

d^k 2 

Here we take into account that Lagrangian is defined within the accuracy 
of the total derivative on time. Then the effective Hamiltonian Tig// = 
^kRk + 'f^kik — ^eff Corresponding to the quark-diquark structure of a baryon 
is given by 

neff = E ^ + f + f + f 4^ + ^^ye - e (ER) - 7 (EO , (102) 

The Hamiltonian for baryons Eq. (102) looks like the one for mesons (44) be- 
cause we consider basically the string between quark and diquark. Therefore 
the calculations of mean coordinates and electromagnetic polarizabilities is 
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the same [31,36]. But in the case of baryons there are more parameters and 
the analysis is more comphcated. 

With the help of Eq. (101) the effective Hamiltonian (102) is rewritten 

as 



H 



eff 



E 



2 

rrij 
2^ 



M 
+ — + 



2M 



n + -(RxH)+7(exH) 



1 

^ 2^ 



7r + ^(exH)+7(RxH) 



+ (76Ve-e(ER)-7(E0, (103) 



where (^ x H)^ = emnkCmHn- The mass of a baryon J^{nj) is defined here 
as a solution to equation 



(104) 



In according to the Noether theorem, the momentum 11 is conserved and we 
can put R = n = in Eq. (103). To find the solution to Eq. (104) we can 
use the substitution TTk — —id/d^k- Then the mass of a baryon is 
given by 

Mil^j) = E ^ + Y + E, H), (105) 



where e(/Xj, E, H) is the eigenvalue of the equation 

-l2 o 



1 

2/1 



■ d A ,^ 

X H) 



+ ^(exH)^ + a6v^-7(E0^$ 



= e(/x„E,H)$. (106) 

The term (7^/(2M)) (^ x H)^ in Eq. (106) is due to the recoil of the string. 
In nonrelativistic models the effect of the recoil was studied in [39] (see also 
[22]). As we neglected the spin of baryons here, there is no interaction of 
spin with the external magnetic field. It is not difficult to take into account 
such interaction. Eq. (106) is like the equation for mesons [36]. 
If E = H = 0, we arrive from Eq. (106) at 



(107) 
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Equation (107) has the same structure as Eq. (51). The eigenvalue of Eq. 
(107) is 



M 



nl/3 



2/i3 (ah + 1^2) 



(108) 



The condition of the minimum of the baryon mass (105) {dAi{fij) / fij = 0) 
at mi = 777.2 (and /^i = 112), with the help of Eq. (105) gives the dynamical 
mass of a diquark: 



a 77 



3/4 



(109) 



This value is different from one [17] obtained for large angular momentum. 
Using Eq. (109) we arrive from Eq. (105) at the expression for the mass of 
a baryon (see also [17]): 



Miiij 



777i + 4m| (0) 



(110) 



To estimate the baryon mass, the value of the string tension cr = 0.15 GeV^ 
win be used [17]. Neglecting the small current masses of quarks m.j wc find 
from Eqs. (109), (110) the mass of a diquark for 77 = 1: nf ^ = 320 McV and 
the nucleon mass : Al(//j) = 1.28 GeV [37]. This value of the nucleon mass 
is a little greater then real nucleon mass because spin-spin and spin-orbit 
forces were omitted. 

Now we consider the mean relative coordinates of nucleons on the basis 
of the virial theorem which gives (as for the case of mesons [36]) the mean 
potential energy (U) = 2(T), where (T) is the mean kinetic energy. Then 
using the relation (T) -|- (U) — e(//j), we arrive at 



(U) = ^6(/.,) = ^(2/.)-V3(6a)^/^a(77 



^2/3, 



(111) 



Comparing Eq. (Ill) with the relation (U) = ha{\f^) gives the following 
expression 



{^e) = -{2^ha)-''''a{n) 




(112) 
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In accordance with Eq. (97) the size of the nucleon is characterized by the 
value I z*^'^) — (l/2)(z'^^) + tP'^) |~| \. Introducing the notation r;, = 6^, 
from Eq. (112) we have 

The same expression was found for mesons (see Eq. (70)). For the quark- 
diquark system the string tension coincides with those of mesons and there- 
fore the quark-diquark system has approximately the same size as mesons. 
Using cr = 0.15 GeV^ and a(l) = 2.2281 we find the mean size of the nucleons 

(y^) = 0.84 fm. (114) 
The experimental value of the charge radii of the proton and neutron are 

^ = QMfm [47], 

(-0.113 ± 0.003) /m [48]. 

In accordance with Eq. (89) the center of mass of the quark-diquark system 
is situated in the center between quark q'^^^ and diquark (c[^^\q^'^^^ and the 

mean radius of a nucleon is (1/2) (y^) — 0.42 fm which is the reasonable 
value. It should be noted that charge radii of hadrons are defined from 
electromagnetic formf actors. 



8 Electric Polarizabilities of Nucleons 

To evaluate electric polarizabilities of nucleons we consider the case when 
H = 0, E 7^ 0. It is possible to assume as an approximation that E || ^ 
[37] , i.e. external electric field is parallel to the string which connects quark 
dth diquark (q^^\ q'^'^^^ ■ So we neglect the rotation of the string. It is 
justified only for the ground state when the orbital quantum number I — 0. 
Introducing the effective string tension 

^e// = ^ - (115) 

we arrive from Eq. (108) at the eigenvalue 

e{i,„E)^{2i,)-'/'{ba,fff'a{n). 
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Prom (109), (110), by neglecting the small terms containing the current 
masses we find the mass of a baryon 



an 



3/4 



(116) 



which depends on the external electric field. Inserting Eq. (115) into Eq. 
(116) and expanding it in a small electric field one yields 



M{i^j,E) 



a(n) 



3/4 



8 



(117) 



where Eq. (109) was used and g = ei + 62 — 63. We write here only terms 
of the expansion 7Vl(//j,E) in small electric field up to E'^ . The first term 
in Eq. (117) gives the mass of a baryon. The second one is connected with 
the potential energy of a dipole moment of quarks in the external electric 
field V — — dE. From Eq. (99) when the center of mass coordinate R = 
and I ^ I ^1 I, the dipole moment of quark-diquark system is d ~ 7^ = 
(7/6)rfe. Comparing the potential energy of a dipole f/ = — (7/6)rfei? (at 
E II r) with the second term of Eq. (117): -{qj ^)\a{n)l'ifl^E we arrive 
at the expression for the mean relative coordinate = (2/y^)[a(n)/3]^/^ 
which coincides with Eq. (113). Here we define the electric dipole moment 
of quark-diquark system more precisely as compared with the letter [37] and 
as a result the mean relative coordinate of a baryon Eq. (113) coincides with 
that for a meson. The third term in Eq. (117) describes the potential energy 
due to the electric polarizability of a baryon (see Eq. (155) in Appendix). 

From Eq. (117) by comparing the quadratic term in E with Eq. (155) 
we arrive at the electric polarizability of a baryon: 



a 



4 



an 



3/4 



r-3/2 



(118) 



Let us consider the estimation of the electric polarizability for the proton 
p = uud. There are two possibilities for a proton as a quark-diquark system: 

a) the quark g*^^) = d and diquark (g^^^g*^^)) = (mm), so the electric charges 
Ci = 62 = (2e)/3, 63 = — e/3 and parameter g = ei + 62 — 63 = (5e)/3; 

b) the quark g^^^ = u, diquark (g^^-'g^^^) = {ud) and the electric charges 
ei = 63 = (2e)/3, 62 = —e/3 and parameter g = ci + 62 — 63 = —e/3. It 
should be noted that there are no permutations of quarks here which occur 
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in the Green function Eq. (79) due to the Pauh principle. Using the value of 
the string tension a = 0.15 GeV^ and a(l) = 2.3381 from Eq. (118) we find 
the static polarizability of a proton in Gaussian units for two cases 

a) ap = 5.56 x 10"^ /m^ b) = 0.22 x 10"^ fm^. (119) 

From Eq. (161) (see Appendix) Aa^ = (4.5 ±0.1) x lO^'^ f'nv'-i and the total 
electric polarizability of a proton becomes (see Eq. (158)) 

a) ftp = 10 X 10"^ /m^ h) = 4.7 x 10"^ fm^. (120) 

The configuration of a proton a), when the quark q^^^ — d and diquark 
(^g(i)g(2)^ = [uu) is more favorable as the experimental values of electric 
polarizability are 

a^^P = (10.9 ± 2.2 ± 1.3) x 10"^ frn^ [49], 
a^^p = (10.6 ± 1.2 ± 1.0) X 10"^ fm^ [50], 

a;^P= (9.8±0.4±1.1) X lO'V^^ [51]. 

So in the case a) we have a good agreement with experimental data. 

For the neutron n = udd there are also two possibilities: a) the quark 
g*^'^^ = u and diquark (g*^^)g*^^)) = (dd), so the electric charges ei = 62 = — e/3, 
63 = (2e)/3 and parameter q — ei-\-e2 — ez — — (4e)/3; b) the quark q^"^^ = d, 
diquark {q'^^^q^'^^) — (ud) and the electric charges 63 — 62 — —e/3, Ci = (2e)/3 
and parameter q — ci + e2 — — (2e)/3. Inserting these parameters into 
Eq. (118) one gives 

a) an = 3.56 x 10"^ /m^ b) an = 0.89 x 10"^ fm^. (121) 

For the neutron, in accordance with Eq. (161) Aq;„ = 0.62 x 10~^ fm^ and 
the generalized electric polarizability of a neutron is given by 

a) an = 4.2 x 10"^ /m^ b) a„ = 1.5 x 10"^ fm^. (122) 

The experimental situation for a neutron is more complicated as there are 
different experimental data: 

a^^P = (0.0 ± 5) X 10-^ fm^ [48], (123) 

a^^P = (12.6 ± 1.5 ± 2.0) x 10"^ fm^ [52]. (124) 
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The recent experimental data (123) are close to both cases but the case b) 
with the quark q'^^^ = d and diquark (q^^^q'^'^^) = {ud) is more favorable. 
Other experimental data [52] arc close to the case a) in Eq. (122) with quark 
q^^^ = u and diquark (g'^^^g'^^)) = {dd) but the magnitude (122) a) is small. 
This experimental discrepancy does not allow us to choose reliably one of 
the possibilities: a) or b). The value of q;„ in the situation a) in Eq. (122) is 
close to one obtained in the oscillator nonrelativistic quark model [53-55,22]. 
The results of calculations in the framework of the dispersion sum rule and 
CHPT Eq. (167) are closer to the case a). 



9 Diamagnetic Polarizabilities of Nucleons 

In according to Eq. (103) for calculating the magnetic polarizabihty of nu- 
cleons one needs to compare it with Eq. (155). The effective Hamiltonian 
(103) (at E = 0, R = n = 0) can be cast into 



Heff — Tio + H 



intt 



Hi 



int 



2/x 1 8// 2M I 



[(exH)]% 



(125) 



where — —iekmn^mdn is the angular momentum and dn —d/d^n- Without 
loss of generality we can choose the direction of the magnetic field on the 
third axis, i.e. H = [0,0, H). Then the Hamiltonian of an interaction of 
quarks with the magnetic field is given by 



(126) 



where L3 = i (^2^1 — Ci92)- Considering the small external magnetic field, 
the perturbative theory can be applied. Using the perturbative method [56], 
within the accuracy of the second order, one arrives at the shift of the energy 



ASr, 



{n 



2^ ^^\8n^2M 



H'^i^ sini? 



n) 
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^ ^ I (nM - (A/f L3) / (2/.) I n) ^^^7) 

where d is the angle between coordinate ^ and magnetic field H. If the first 
term in Eq. (127) is not equal to zero then the second and third terms are 
smaller, and the main contribution to the energy comes from the first term. 
This occurs when the orbital quantum number Z > 0. In the case I — 0, the 
shift of the energy due to the interaction with the magnetic field is defined 
by the second term in Eq. (127). 

After averaging in Eq. (127) and taking into account the equation 



sin^ = \ 



we find for the ground state when I — 0, the shift of energy 

where {^'^) = (0 | | 0), and | 0) means the wave function of the ground 
s-state. We took into account that the first and the third terms in Eq. (127) 
equal zero because L3 | 0) = 0. Here we ignore the spin interactions of 
a baryon with the external magnetic field and therefore only diamagnetic 
polarizability can be defined from Eq. (128). Using the definition of the 
relative coordinate r;, = and comparing Eq. (128) with Eq. (155) gives 
the diamagnetic polarizability of a baryon 

^'-''--i^r + t-MrM- (129) 



As required, the diamagnetic polarizability is negative and Eq. (129) is like 
the Langevin formula for the magnetic susceptibility of atoms. The similar 
expression was derived in [36] for mesons. Taking into account Eqs. (92), 
(109), expression (129) is rewritten as 

'3*° = -(t + '')||. (130) 

where M = jJ^i'^ + 1^2'^ + fJ^^'^ = '^^jJ^f^ = 2^/ct [a(n)/3]^^'^. To calculate the value 
of for nucleons we can use the theoretical magnitude of the mean-squared 
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relative coordinate (r^) or experimental data for the size of a nucleon. The 
first way is preferable. Taking into account Eq. (113), the value of M and 

using the approximate relation ( r^) ~ (Wr^)^, Eq. (130) transforms into 



pdta 



12(t3/2 



a{n) 



(131) 



For a proton with the more favorable configuration a), when the quark q^^^ = 
d and diquark (q^^^q^"^^) = (uu), parameter q = (5e)/3 and Eq. (131) (at 
a — 0.15 GeV^) takes the value 

pdia = _8 X 10-4 fm\ (132) 

This quantity is greater than that found in the nonrelativistic quark model 

[22] . To calculate fjP"'^"' in the present approach one needs to take into account 
the interaction of spins of quarks with the magnetic field. For estimation of 
the total magnetic polarizability of a proton we use the contribution (3^^"' 
and in accordance with Eqs. (159), (162) we find 

Pp^{5±3)x IQ-^ fm^. (133) 

This quantity is close (within two standard deviations) to the experimental 
value [52] which was extracted from the measurements of the cross sections 
of the Compton scattering on hydrogen: 

= (2.9 T 0.7 T 0.8) X 10-^ fm^. (134) 

The value (133) is also in agreement with quantity (167) found in other 
schemes. The value Pp is small due to the partial cancellation of the posi- 
tive paramagnetic polarizability /?^"''" with the negative diamagnetic polar- 
izability Pp^"'. The theoretical evaluation, interpretation and experimental 
extraction of the total magnetic polarizability /3p is difficult because it is 
small. That is why this quantity is model dependent in different schemes. 
The approach considered allows us to improve the accuracy by using the 
perturbation in the spin interaction. The next step is to take into account 
such corrections. 

For a neutron with the configuration a) where the quark q^^^ = u and 
diquark (g^^^g*^^)) = (dd), the parameter q = — (4e)/3 and Eq. (131) gives 

l^dia ^ _5_4 X 10-4 y^3_ ^^35^ 
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Using the paramagnetic polarizability Z?^"*^" = (13 ± 3) x 10~^ fm^ [57] for 
a neutron (see Eq. (162)) we get from Eq. (135) the generalized magnetic 
polarizability of a neutron 

= (7.6 ± 3) X 10-^ fm^ (136) 

which is in agreement within two standard deviations with the experimental 
quantity [52] 

= (3.2 T 1-5 T 2.0) X 10"^ fm\ 

According to Eqs. (120), (122), (133), (136) the sum of nucleon polarizabil- 
ities in our approach are given by 

ap + Pp^ (15 ± 3) X 10-^ /m^ 

an + (3n = (11.8 ±3) X 10-^ fm^. (137) 

Values (137) arc close to the reliable results found on the basis sum rule 
(166). As pion degrees of freedom play a very important role in the nucleon 
electric polarizability, one needs to describe pions in this scheme. In the 
approach considered here the accuracy of the values of electromagnetic po- 
larizabilities can be improved by taking into consideration spin interactions 
as perturbations in accordance with Eq. (87). 

10 Cluster Expansion 

In this paragraph we go into Euclidean formalism and follow closely the 
works [16,17,61,66]. Some important characteristics of strong interactions 
depend on the Wilson integral averaged with certain weights. So, spin-spin 
and spin-orbit quark interactions can be expressed through the Wilson loop. 
It is difficult to calculate the Wilson loop integral in the general case of 
arbitrary contours due to the nonlinear character of gluonic fields. However, 
it is possible to use the powerful method of cluster expansion [64] . To use this 
method it is necessary to accept the stochastic nature of the QCD vacuum 
[16]. Lattice data confirm this conception [29, 61] and give the correlation 
length Tg ~ 0.2 fm. 

The method of field correlators allows us to evaluate approximately the 
Wilson loop using a cumulant expansion. MFC takes into account both 
perturbative and nonperturbative interactions of gluons and therefore it is 
applicable for small and large distances. 
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The Wilson integral is expressed via gauge-invariant and Lorenz covariant 
correlators of strengths of gluonic fields. To make this expansion one needs 
to use the non-Abelian Stokes theorem (25). 

Let us introduce the generating functional for field correlators of gluonic 
fields (see [65]) as follows 

Z[J] = (exp (^g J dxJi^i,{x)Fi^^,{x,xo)^), (138) 



where Ff^i,{x,XQ) = ^{xo,x)Ff^^{x)^{x,xo). Using the Taylor expansion in 
J^i/ we arrive at 

Z[J] = 1 + f dXn - ■ ■ f dXnJ ^LivA^l) ■ ■ ■ J^iuVui^n) 

n=l ^- J 

X {FiiiuA^l^Xo) ■ ■ ■ F^r.'^^{Xn,Xo)), (139) 

so that the field correlators (averaged strength tensors over gluonic fields) 
are given by 

(io)~'^S^Z\J] 

{F^,.A^uXo) ■ ■ ■ F^^^^xn^xo)) = y \j=o . (140) 

It is convenient instead of expressions (140) to use irreducible correlators - 
cumulants [64] in accordance with relationship 

{{F^Mx,,xo) ■ ■ ■ F^„,„{xn,xo))) = XT l^=° • 

Prom Eq. (141) we come to the expansion of the generating functional Z[J] 
via cumulants: 

Z[J] = exp[l + f dxi - ■■ f dXnJ,j.^v^{xi) ■ ■ ■ J^^vS^n) 

^ n=l f^- J J 

X {{F^l^uAxl,XQ) ■ ■ ■ F^„^„{xn,xo)))y (142) 

Making variations of Eqs. (139) and (142) and setting J^,^ = one arrives at 
the first three relations [64] (we omit here argument xo implaying F^,/(x, xq) — 
FU^)): 

{{F,^{x))) = {F,^{x)), 
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ix2)F^,,,,{x3))) = {F^,^^{xi)F^^^^{x2)Ff,^^^{x3)) 

~{F^iuAxi)F^2y2i.X2)){F^;vAx^)) - {Fi,^uAxi)){F^2u^{x2)F^^^.^{x^)) 

-{F„^yAxi)F„zVi{xz)){F^2V2{x2)) + {F^rvr{xi)){F^^y^{x2)){F^^y:,{x^)) 

+ {F,,,A^,)){F,,,,{x^)){F,,,,{x2)). (143) 

It should be noted that there are two types of ordering in cumulant expan- 
sions [17], [64]. According to the definition of [17], there are no terms in Eq. 
(143) hke 

{F,,^uA^i)F,,^vz{xz)){F^2U2{x2)), {Fi,^uAxi)){F^^^^{x3)){Fi,^^^{x2)), 

which violate path-ordering prescription. For the Gaussian approximation 
when we retain only the bilocal cumulant both definitions coincide. The 
average value (F^,y(a;)) should vanish, because in the stochastic vacuum there 
is no the definite direction. Otherwise the Lorentz invariance will be broken 
by vacuum. Applying the cumulant expansion (142) to the average Wilson 
integral (24) one arrives at 

{W{C)) = ±trPe^p(f: / da,,^,{x,) ■■■ f da,^M 

X {{F^,,vAxl,Xo) ■ ■ ■ F^„^„{xn,Xo)))^, (144) 

where P an ordering operator on the coordinate x of the matrix Ff^:y{x,xo), 
and we imply averaging on gluonic fields with the standard weight (see (5)). 

The expansion (144) makes sense when it converges. As we assume that 
Eq. (144) is valid for any g (even for large distances), cumulant expansion is 
nonperturbative. Lattice data showed [29, 66] that the correlators decrease 
when the distance between two points increases and the expansion (144) is 
justified. It is known that in a Gaussian random process, the higher cumu- 
lants vanish, and only the quadratic cumulant does not equal zero. With 
good accuracy the QCD vacuum can be considered as a stochastic ensemble 
of gluonic fields [66] . This means that the main contribution to the Wilson 
loop expansion (144) comes from the bilocal cumulant 

{{F^Ax,Xo)F^p{y,Xo))) = {{F^,{x)^x,y)F^p{y)^y,x))). (145) 



41 



Higher cumulants in Eq. (144) give small corrections with {{F^i,Y)Tg <S 1 
[16, 17], and can be neglected (Gaussian approximation). The Gaussian dom- 
inance was supported by lattice calculations and the hadron phenomenology. 
Thus we come to the Gaussian stochastic model of the QCD vacuum [29, 
66]. The fundamental input in such a model is the Gaussian correlator. But 
there is a problem; what kind of field configurations (classical or quantum) 
play the dominant role to form the definite behavior of field correlators. 

However, there will be the dependence of Wilson loop integral (144) on 
the point xq and the shape of the surface S. We recall that the point xq and 
contour C belong to the surface E. The total expression (144) with higher 
cumulants does not depend on the point Xq and contour C. To reduce the 
dependence of a cumulant on the E, one should choose the surface with the 
minimal area [17]. Then the area law of the Wilson integral with large loops 
(the size i? ~ 1 fm) takes place [17]. To get rid on the dependence of the 
Wilson loop on xq one should have the condition [17] 

I x'^^ - x^^^ ]<| x^^ -xo\, I x^^^ -xo\ . 

Thus, at small correlation length Tg <^| x*^*) — xq \ one can neglect the de- 
pendence on Xq. In the Gaussian approximation the Wilson integral (144) 
becomes 

X {{F^,{x,xo)F^p{y,xo)))^, (146) 

where we imply that x — x(^), y — y{C). Thus, for a size of the Wilson 
loop Tg the bilocal cumulant (145) entering Eq. (146) does not depend 
on Xq. As a result it can be expressed via two scalar functions Di and D2 
[16,17]: 

r 

— {{Fi,^{x)^{x,y)F^p{y)^{y,x))) = iNcii^na^up - ^i,o5ya)F> {x - y) 



1 

+2 



d , , 

X - y)aSul3 -{X - y)j35r, 



dXf, 



+ -^{{x-y)(i5na-{x-y)JnP^ -Di(a;-7/)|, (147) 
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where I^c is a unit matrix Nc x Nq-, and functions D and Di are invariants 
of the renormahzation group. As the cumulant (147) is proportional to the 
unit matrix Inc) ^he ordering operator P can be omitted in Eq. (146). 
Substituting Eq. (147) into Eq. (146) one arrives at Eg. (29) with the string 
tension [17]: 

a j (fxD{x). (148) 

So, a is an integral characteristic of the Gaussian correlator. Only function D 
enters Eq. (148) defining the string tension but both functions D and Di give 
the contribution to the perimeter term ■jL{C) in Eq. (31) [17]. Higher cumu- 
lants add small corrections to the string tension (148). The functions D, Di 
in the nonperturbative regime were derived in lattice experiments [29,66] and 
are given by 

D{x) = D{0)exp(-^-^), Di(a;) = L'i(0)expf-^y (149) 



\ J.g / \ ±g / 

where L'(O) = 0.073 GeV^, Di(0) = (1/3)D(0), Tg = 0.2 0.3 fm. Values 
D{0), Di{0) can be connected with the gluonic condensate [4] 

^((f;,(0))V 0.012 Gel-^ (150) 

so that D{0) ^ 2as((i^;?^(0))Y Thus, MFC is a development of the QCD 
sum rule method. In the perturbative regime at small distances, the contribu- 
tion from gluonic exchange (giving the Coulomb potential) to the D{x), Di{x) 
is ofO{g') [17] 

16q;s 
3nx^ 

The colour Coulomb interaction is dominant below 0.3 fm. This region is 
important for bottomonium (hb state) and charmonium (cc state) where the 
sizes of ground states are around 0.2 fm and 0.4 fm, respectively. For cc 
and ss systems both perturbative and nonperturbative interactions should 
be taken into account, and only two D, Di functions can be used. It should 
be noticed that nonperturbative sector influences the perturbative part and 
vice versa. So, two functions, D, Di can be represented as [17,66] 

D{x) = L>^^'-*(x) + L>^^(x), D^{x) = D^^'\x) + D^''{x). (152) 



D{x) = Q, D,{x) = ^^. (151) 
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At — > the perturbative part of the function D has the form D^^^^{x) ~ 
q; > 0. The nonperturbative part is normahzed by the gluon condensate 



D^^(O) + Df^(0) = (l/24iVc)((F;,(0))') = (7rV6iVc)(0.14±0.02) GeV\ 



The value of the gluon condensate was taken from the lattice data. 

The strong nonperturbative chromoelectromagnetic fields of the QCD 
vacuum cause nonperturbative shift of the energy density of the vacuum due 
to the scale anomaly [4] : 



We know that the Gell-Mann - Low function, [3{as) < due to asymptotic 
freedom at small a^. So, the nonperturbative shift of the vacuum energy 
density is energetically favorable. 

Because only function D[x) contributes to the string tension (148), it is 
responsible for a condensation of monopoles in accordance with the 't Hooft 
and Mandelstam scenario of confinement. The lattice calculations gave the 
width of the string / = 0.4^0.5 /m, so that the correlation length Tg ~ (1/2)1. 
Thus, the radius of the string is close to Tg. Therefore another interpretation 
of Tg is the thickness of the confining string. In the case — >^ we have 
pure string limit. The reader can find more about MFC in [17,66]. 

11 Conclusion 

The QCD string theory allows us to obtain the effective Lagrangians for 
mesons and baryons and estimate the mean squared radii and electromag- 
netic polarizabilities of hadrons. These quantities were derived as functions 
of the string tension which is a fundamental variable in this approach. It is 
not difficult to calculate the electromagnetic polarizabilities of excited states 
of hadrons using this approach. For that we should take the quantum num- 
bers = 1, / = (n = 2) and evaluate the mean relative coordinate 
in accordance with Eq. (113). Then Eqs. (118), (130) give the required 
polarizabilities. To have more precise values of the hadron electromagnetic 
characteristics one needs to take into account spin corrections. Especially it 
is important for light pseudoscalar mesons (tt, K mesons). In principle it is 
possible to receive the spin-orbit and spin-spin interactions using the general 




16Q!s 



(153) 
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expressions (23), (83). The model of baryons as a quark-diquark system is 
very similar to the approach for mesons. The mean size and electromag- 
netic polarizabilities of a proton and neutron are in reasonable agreement 
with the experimental data. The more favorable combination for a diquark 
is (uu) for a proton and {dd) for a neutron. In this case theoretical values 
for electromagnetic polarizabilities are close to experimental data. 

It should be noted that the quark-diquark ansatz used here generates an 
electric dipole moment, and a result, there is a contribution to the energy 
linear in the electric field which is like the linear Stark effect. This requires 
the deeper understanding the properties of a nucleon under the symmetry 
transformations (for example parity). 

Spin interactions of quarks treated as a perturbation were neglected here 
but we took them into account by using the paramagnetic polarizability of 
a nucleon due to the A contribution. As an approximation it is justified 
because in Isgur-Karl model of baryons [58] spin-orbit splitting are much 
smaller than expected from one-gluon-exchange matrix elements (spin forces 
were also discussed in [59]). It was also shown [60,61] that the contributions 
from the Coulomb and spin-spin interactions cancel each other. Nevertheless 
all parameters should be taken from the approach and then compared with 
empirical data. Therefore one needs to calculate the paramagnetic polariz- 
ability of a nucleon in our approach. 

Implying small spin-orbit forces in baryons we come to K* and quark- 
diquark system correspondence (see [36]). The present approach can be ap- 
plied for studying any baryon (see also [60]). 

The Nambu-Jona-Lasinio (NJL) model [13,14] having a basis in the frame- 
work of QCD [12] describes the chiral symmetry breaking but not the con- 
finement of quarks [17,15]. Besides this model has free parameters and the 
calculated polarizabilities of mesons [62] are parameter dependent. 

The instanton vacuum theory developed in [11] does not give the con- 
finement of quarks phenomenon [17]. This theory is like the NJL model [15] 
and takes into account only the chiral symmetry breaking. Therefore the 
calculation of the meson electromagnetic polarizabilities on the basis of the 
IVT gave the similar result [63] as in the NJL model. 

The electromagnetic polarizabilities of nucleons found are close to the 
values calculated in the framework of the dispersion sum rule [71] and the 
chiral perturbation theory (CHPT) [72]. 

To improve the accuracy of calculations of electromagnetic polarizabilities 
of nucleons one should take into account pion degrees of freedom because 
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the pion cloud contributes substantially to electromagnetic properties (see 
Appendix). 

All this shows that the theoretical evaluation of the charge radii and the 
magnetic polarizabilities of hadrons is possible on the basis of a good descrip- 
tion the chiral symmetry breaking and confinement of quarks in the frame- 
work of the QCD string theory but with some approximations and model 
assumptions. Naturally this theory was derived using the non-perturbative 
QCD, i.e. first principals of QCD. 



In this Appendix we define nucleon electromagnetic polarizabilities in dif- 
ferent approaches. The static electric a and magnetic (5 polarizabilities 
of a hadron characterize its internal structure. These quantities are low- 
energy characteristics of a hadron and therefore they are defined in the non- 
perturbative regime of QCD. In external electromagnetic fields the hadron is 
deformed and induces dipole electric and magnetic moments 



The acquired dipole moments (154) of a hadron give the contribution to the 
potential energy as follows 



We use Gaussian units here. The cross section of scattering of a photon 
on hadrons (Compton scattering) depends on the electromagnetic polariz- 
abilities due to Eq. (155). Thus, low-energy Compton scattering can give 
important information about the internal structure of hadrons. Below we 
discuss the nucleon electromagnetic polarizabilities. The Compton ampli- 
tude for scattering a photon on a nucleon at low energies in the nucleon rest 
frame (the laboratory system) is given by [22] 



/(77V ^ -iN) ^fB + ujuj'a{ee') + p{k x e)(k' x e') + ^(cu^), (156) 



where (k, io;), e and (k', io;'), e' are four momenta and polarization vectors 
of the incoming and outgoing photon, respectively. The electromagnetic po- 
larizabilities a, P occurring Eq. (156) are so-called Compton polarizabilities. 



12 Appendix 



D = aE, M = (3U. 



(154) 



(155) 
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Below we will discuss the interrelation between static and Compton polariz- 
abilities. The fs is the Born amplitude which is due to the nucleon electric 
charge eZ and anomalous magnetic moment k of a nucleon. The Born am- 
plitude Jb includes the Thomson amplitude /r which is energy-independent 
and satisfies the well-known low-energy theorem: 

where Mn is the nucleon mass. The Thomson scattering amplitude (157) 
corresponds to a pointlike particle and the total amplitude (156) includes 
structure parameters - electromagnetic polarizabilities a, /3. 

In the nonrelativistic approximation the generalized electric and magnetic 
polarizabilities which are extracted from measurements of the Compton scat- 
tering cross sections are given by [22] 

a — a + Aa, 

„ . , 5: U^iiaHLf. (158) 

/^- = 2i:^!iiMli^. (159) 

where | n) is the excited state of a hadron, and are the third projec- 
tions of the electric and magnetic dipole operators, respectively. The excited 
state I n) includes the meson-baryon intermediate states. In nonrelativistic 
calculations the Aa is due to form factors and depends on the electric charge 
radius. Eqs. (158), (159) are similar to the usual quantum-mechanical for- 
mulas. Eq. (158) was also confirmed by the quasiclassical calculations of 
pion polarizability in the framework of the instanton vacuum theory [63]. 
In the completely relativistic approach the physical meaning of a is more 
complicated [22] . In [22] the correction Aa is given by 

where is the electric radius and the magnetic moment of the hadron is 
given hy fjL— (ck) / (2Mjv) . 
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For the proton and neutron, using the experimental values of electric radii 
and magnetic moments, the term (160) becomes [22] 

Aap = (4.5 ± 0.1) X 10-^ /m^ Aq;„ = 0.62 x 10"^ fm^ (161) 

The strong magnetic A^A transition gives rise to a positive part of the mag- 
netic polarizabilities. Therefore the main contribution to the paramagnetic 
polarizabihty of a nucleon is due to A(1232) excitation [57] and is given by 

The intermediate state A is an isovector excitation and contributes to the 
proton and neutron so that — P'^. The interesting feature of the experi- 
mental data is that for the proton and the neutron the approximate equality 
(ip + Pp — ctn + Pn, and inequality dp,dn 3> Pp,Pn are valid. The last in- 
equality means that both a proton and a neutron behave like electric dipoles. 
In addition, the large positive paramagnetic polarizabihty of nucleons should 
be canceled by the diamagnetic contribution. The sum of electromagnetic 
polarizabilities can be connected with the total photoabsorption cross sec- 
tion [67,68,22]. Indeed, from Eq. (156) the spin-independent amplitude of 
the Compton scattering at low energies and zero angle and e = e' is given by 

f,^y{u, 0) = + (d + P)u;' + 0(u;'). (163) 

Using the optical theorem [22] 

Imf^^Y{uj,0) ^ ^atotiuj), (164) 

and once-subtracted forward dispersion relation for amplitude (163), one 
finds the relation [67,68] 

1 duj 

{a + p)^— —atotico), (165) 

where a; = (s — Mfj) / {2Mm) is the photon energy in the laboratory sys- 
tem and uithr = ^7r(l + ^7r/(2MAr)) is the pion production threshold (m^r = 
139.57 MeV is the pion mass). Eq. (165) occurs also for spin-0 particles. 
The strict equality (165) allows us to extract the sum of electromagnetic po- 
larizabilities from the experimentally measured total photoabsorption cross 
section. The result is [69,70] 

{dp + pp) = (14.3 ± 0.5) X 10"^ /m^ 
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(a„ + Pn) = (15.8 ± 0.5) X 10-^ fm^. (166) 

The most difficult task is to find electromagnetic polarizabilities separately 

(sec [22]). 

It follows from Eq. (165) that at energies u ~ the sum of polariz- 
abilities behave as 0(m~^) because the total photoabsorbtion cross section 
atoti'-^) is finite at a; ~ m^r — > [71]. In the chiral limit, —>■ 0, both po- 
larizabilities a, P diverge as 0{m~^). This was confirmed in the framework 
of the chiral perturbation theory (CHPT) [72]. The same behavior of pion 
polarizabilities was found on the basis of quasiclassical calculations [63] in 
the framework of the instanton vacuum theory. Thus, the leading terms of 
electromagnetic polarizabilities are determined by CSB. 

The contribution of low-lying intermediate states, where soft pions play 
an important role, into the integral of Eq. (165) was estimated in [71]: 

2 

5 (eg^NV^\ io«,, -,n-4^ 3 
« — — — — 77 — — 13.6 X 10 jm , 



12m^ V 871 



N 



^ ^ J_ f 1.4 X 10-^ /m^ (167) 

where the strong pion-nucleon coupling constant Qt^n = 13.4 is connected 
with the charged pion decay constant = 93 MeV by the Goldberger- 
Treiman relation g^j^/M^- = gA/fn, where gA is the axial- vector coupling 
constant. The values (167) found for the leading contribution in the chiral 
limit, and therefore a, P are the same for the proton and neutron. The 
next corrections 0{uj/Mm) to polarizabilities, due to the nucleon recoil, are 
different for a proton and neutron [71]. The result (167) is in agreement with 
CHPT calculations [72]. 

The nucleon polarizabilities were also evaluated on the basis of nucleon 
soliton models in [73-79] . 
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